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ABSTRACT 



The TRACER instrument ( "Transition Radiation Array for Cosmic Energetic 
Radiation" ) has been developed for direct measurements of the heavier primary 
cosmic-ray nuclei at high energies. The instrument had a successful long-duration 
balloon flight in Antarctica in 2003. The detector system and measurement 
process are described, details of the data analysis are discussed, and the individual 
energy spectra of the elements O, Ne, Mg, Si, S, Ar, Ca, and Fe (nuclear charge 
Z=8 to 26) are presented. The large geometric factor of TRACER and the use 
of a transition radiation detector make it possible to determine the spectra up 
to energies in excess of 10 14 eV per particle. A power-law fit to the individual 
energy spectra above 20 GeV per amu exhibits nearly the same spectral index 
(~ 2.65 ± 0.05) for all elements, without noticeable dependence on the elemental 
charge Z. 
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Introduction 



The energies of cosmic rays observed near earth extend over a very wide range, from 
below 10 8 eV to more than 10 20 eV per particle. Up to about 10 10 eV per amu, the particle 
energies and intensities are significantly affected by solar modulation, but for the nearly 
ten remaining decades, the arriving particles are believed to represent the ambient cosmic- 
ray population in the local galaxy. Over this range, the total cosmic-ray i ntensity (the 



all p article" differential energy spectrum, see, for instance, the compilation of ICronin et al. 



( 119971 )) decreases monotonically by more than 30 orders of magnitude, initially as a power law 
oc E~ 2 - 7 , steepening slightly to E~ 3 at the "knee" above 10 15 eV per particle, and exhibiting 
small changes in slope again (a "second knee" and an "ankle") in the 10 18 - 10 19 eV region. 
Besides these relatively minor changes, the overall spectrum is remarkably featureless, even 
though a variety of processes may contribute to the cosmic-ray flux. The general consensus 
is that the particles below the knee are generated inside our Galaxy and are contained by 
the galactic magnetic fields for millions of years, while at the highest energies (> 10 18 eV) 
the gyration radii become so large that galactic containment is no longer effective, and that 
the cosmic rays may then be of extragalactic origin. 

The all-particle energy spectrum refers to the overall spectrum of cosmic rays without 
differentiating the individual components. To obtain deeper insight, details of the composi- 
tion of the cosmic-ray particles must be studied. For instance, very accurate measurements of 
not only the elemental, but also the isotopic composition are now available at low energies (< 
10 9 - 10 10 eV per particle). These have determined such important quantities as the average 
containment time of cosmic rays in the galaxy fro m observations of rad i oactiv e clock-nuclei 
(r ~ 1.5 x 10 7 y at these ener gies (first reported by lGarcia Mufioz et al.l (119751 ); more recent 
work by lYanasak et al.l (120011 )). and they have excluded fresh super nova e.jecta as source ma - 
terial of cosmic rays from the abundances of Co- and Ni-isotopes ( jWiedenbeck et al.lll9991 ). 
At relativistic energies only the abundances o f the elemental species, w ithout isotopic detail, 
have been accessible to measurements (e.g. lEngelmann et al.l (119901 ) ) . but even these en- 
counter increasingly severe systema tic and statis t ical u ncertainties as the energy increases. 
Very few direct observations (e.g. iMiiller et al.l (Il99ll )) have provided composition detail 



with single-element resolution above 10 13 eV per particle. The work described in this paper 
represents an attempt to improve this situation, and to determine the energy spectra of 
individual elemental species up to much higher energies. 

The current paradigm of the origin of galactic cosmic rays postulates first-order Fermi 
acceleration in interstellar shock fronts from supernova (SN) explosions as the most p otent 
contri butor to the cosmic-ray flux below the knee. This mechanism, first proposed by iBell 
(119781 ) . must be very efficient indeed, since as much as ~10% of the total kinetic energy 
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released in supernova explosions is required to sustain the cosmic-ray flux. The SN shock 
acceleration mechanism predicts a source energy spectrum in the form of a power law cx E~ r , 
with a spectral index T « 2.0 for strong shocks. Such a spectrum is much harder than the 
E~ 2 - 7 spectrum that seems to be typical for most primary cosmic-ray nuclei below the knee. 
This difference in spectral slope can be explained if the propagation and containment of 
cosmic rays in the galaxy are energy- dependent. Such a behavior had indeed been inferred, 
even before the SN shock acceleration model was formulated, on the basis of measurements 
of the "L/M ratio", i.e., of the abundance of secondary, spallation-produced cosmic rays 
(such as the light elem ents Li, Be, B) relative to their primary parents (such as C and O), 
(IJuliusson et al.lll972l ). One concludes from such measurements that the average amount of 
interstellar matter encountered by cosmic-ray nuclei during their propagation through the 
galaxy decreases with energy. The effect is often parameterized by a propagation pathlengt h 
A, which depends on energy as A oc E~ os for relativistic nuclei (jEngelmann et al.lll990l ). 



The primary cosmic-ray spectrum at the source would then be approximately oc E~ or 
E~ 2 - 2 , close to what the shock acceleration model predicts. This fact provides strong but 
indirect evidence for the validity of the shock acceleration model. 

However, it must be kept in mind that current measurements of the L/M ratio do not 
extend beyond 10 11 eV per amu (or about 10 12 eV per particle). Hence, there is an extended 
region of energies (10 12 - 10 15 eV per particle) where the SN shock acceleration model is 
assumed to be valid, albeit without the benefit of much observational support. The SN 
shock acceleration process is expected to become ineffic ient at energies around Zx 10 14 eV 
(where Z is the atomic number of a cosmic-ray nucleus) (ILaGage &: Cesarskylll983l ). but this 
limit is not observationally confirmed. One might refer to the steepening of the all-particle 
spectrum at the "knee" as evidence, but even then the origin of cosmic rays with energies 
far beyond the knee remains a mystery. 

The SN-shock acceleration model has recently found strong o bservational support th rough 
the detection of TeV gamma rays from shell-type SN-remnants (lAharonian et al.l 12004 ). al- 
though there remains a debate whether the parent parti cles of these gamma rays are electrons 
or, indeed, protons and nuclei (IBerezhko & Volkl|2006l ). 



Clearly, the present observational evidence about high-energy cosmic rays is inadequate 
to provide answers to some of the most fundamental questions about their origin and galactic 
propagation, and there exists an undeniable need for improved measurements at higher 
energies. The TRACER program has been developed to provide direct measurements of the 
energy spectra of the heavier cosmic-ray nuclei through balloon flights above the atmosphere, 
up to energies between 10 14 and 10 15 eV per particle [The acronym TRACER stands for 
"Transition Radiation Array for Cosmic Energetic Radiation"]. Higher energies are currently 
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studied through indirect observations with ground-based air-shower installations. However, 
attempts to obtain composition details from air-shower observations are affected by a number 
of systematic uncertainties. One of the goals of TRACER is to approach the energy region of 
air shower measurements, with the hope of providing some cross-calibration with the indirect 
measurements. 



2. Observational Technique 

2.1. General Principles 

A successful cosmic-ray measurement above the atmosphere must determine at least two 
quantities for each particle: the charge Z and the energy E. Measurement of Z is commonly 
accomplished by utilizing the fact that all electro-magnetic interactions scale with Z 2 . Hence, 
devices that measure the ionization energy loss (dE/dx) such as plastic scintillators or gas 
proportional counters, or that measure the intensity of Cherenkov light, are well suited and 
are used in TRACER. 

The energy determination is a greater challenge, mainly because of the large detector 
areas that are needed at high energy. Extrapolations from lower energies indicate that 
a successful measurement of nuclei heavier than helium requires detectors with exposure 
factors of at least several 100 m 2 sr days in order to approach the 10 15 eV per particle 
region. Such exposures can be accomplished with long-duration (i.e., several weeks) balloon 
flights if the detectors have sensitive areas of a few square meters. A classic calorimeter 
that measures the total energy absorbed in matter after a particle interacts would become 
much too massive if it were to have an area of this magnitude. Hence, TRACER uses 
electromagnetic interactions not just to determine Z, but also to measure the energy (or 
more exactly, the Lorentz-factor 7 Pd E/mc 2 ): measurement of Cherenkov light produced in 
an acrylic Cherenkov counter at low energies (up to a few 10 9 eV per amu), measurement 
of the specific ionization in gases at intermediate energies (10 10 to 10 12 eV per amu), and 
measurement of transition radiation at the highest energies (> 5 x 10 11 eV per amu). 

Figure [T] shows the typical energy responses for the detectors used in TRACER. The 
light intensity of an acrylic Cherenkov counter (refractive index = 1.49) increases quickly 
with energy above a threshold of 0.325 GeV per amu (7 ~ 1.35) until it approaches saturation 
above 20 GeV per amu (7 ~ 22). The ionization loss in gas decreases with energy, reaches a 
minimum at 7 = 3.96, and then increases logarithmically with E. Transition radiation (TR) 
x-rays are generated in a multi-surface radiator (mats of plastic fibers in TRACER) and are 
detected in a gas proportional counter. Hence, the TR detector measures an ionization loss 
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signal, and at sufficiently high energies (> 5 x 10 11 eV per amu), a superimposed signal 
due to the absorption of TR photons. The TR signal rises rapidly with the Lorentz factor 7 
and reaches saturation in the 7 = 10 4 - 10 5 region. The use of TR to me asure the highest 



cosm ic-ray energies was first implemented for the CRN detector in 1985 (IL'Heureux et al. 



19901 ). but remains unconventional. It permits the construction of detectors of very large 
area. As the TR intensity increases with Z 2 , this technique is well suited for measurements 
of the heavier cosmic rays (Z > 3), but not applicable to protons or He-nuclei where the 
signals are affected by large statistical fluctuations. 

Note that the response curves shown in figure [1] for the detectors of specific ionization 
and TR are double-valued: large signals can either be produced by high-energy particles, 
or by particles below minimum ionization. The TRACER instrument must be launched at 
latitudes with low geomagnetic cutoff energies where a significant cosmic-ray intensity arrives 
with sub-relativistic energies. The inclusion of a plastic Cherenkov counter which identifies 
the low-energy particles is then essential in order to remove the degeneracy of ionization and 
TR response. Overall, the counter combination chosen for TRACER measures the energies 
of cosmic-ray nuclei over a very wide range, in excess of four orders of magnitude. 



2.2. The TRACER Instrument 

A schematic cross-section of the TRACER instrument is shown in figure [2j The detector 
elements are rather large in area (206 cm x 206 cm), and are vertically separated by 120 
cm. This leads to an overall geometric fac tor of about 5 m 2 sr. The instrument has been 



previously described (IGahbauer et al.ll2004j ). and a detailed technical description of the de- 
tector elements is currently being prepared for a separate publication. Here, we will present 
a brief summary of the relevant properties of the instrument. 

As figure [2] shows, the individual components of TRACER are, from top to bottom: 

1. A plastic scintillator sheet (BICRON-408 material), with active area 200 cm x 200 
cm, and 0.5 cm thiclfl The scintillator is viewed by 24 photomultiplier tubes (PMT) via 
wavelength shifter bars. 

2. Four double layers of single- wire proportional tubes, oriented in two orthogonal 
directions. Each tube is 200 cm long and 2 cm in diameter. The tube walls are wound from 
aluminized mylar and are about 127 /zm thick. 



1 Due to the interspersed wavelength shifter-bars, the geometric area of these counters is slightly larger 
than the active area. 
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3. Another four double layers of proportional tubes of identical design, but each double 
layer located below blankets of plastic fiber material, which form a radiator to generate 
transition radiation. 

All tubes are filled with a mixture of xenon and methane (at equal parts by volume) 
and are operated in flight at a pressure of 0.5 atmospheres. 

4. A second scintillator of identical design, but read out by just 12 PMTs, and located 
below the lowest layer of proportional tubes. 

5. A Cherenkov counter, with active area 200 cm x 200 cm and 1.27 cm thick, at 
the bottom of the detector stack. The Cherenkov material consists of an acrylic plastic 
(refractive index = 1.49), doped with wavelength shifter material. Signals are read out via 
wavelength-shifter bars connected to 24 PMTs (as for the scintillator on top). 

The two scintillators serve as triggers for the instrument and, together with the Cherenkov 
counter, determine the elemental charge Z of individual cosmic-ray nuclei traversing the in- 
strument. The proportional counter arrays measure the energy (or Lorentz-factor rs E/mc 2 ) 
of the cosmic-ray particles at high energies, and they determine the trajectory of each particle 
traversing the instrument in two projections. The signals in the top four double layers simply 
are a measure for the ionization energy loss idE/dx) of the particle. The small logarithmic 
increase above the minimum-ionization level provides an estimate of the particle energy up 
to Lorentz factors of a few hundred. The lower four double layers of tubes form a transition 
radiation detector (TRD). As the response curves in figure [T] show, the TRD signals are 
identical to the dE/dx signals of the tubes above, except at the highest energies (7 > 400), 
where the rapid increase of the signal with increasing 7 permits an energy measurement. 
The signals of th e proportional t ubes (a total of about 1600) are processed through VLSI 



chips ( "Amplex" . iBeuville et al.l (J1990J)). The finite dynamic range of these devices restricts 
the range of elemental species for which data could be recorded in this flight, from oxygen 
(Z = 8) to iron (Z = 26). 

The TRACER instrument derives its heritage fro m the CRN detector t hat was devel- 



oped by our laboratory in the 1980s for space flight (IL'Heureux et al.lll990l ). In fact, the 
radiator material used for the TRACER-TRD is identical to that of CRN. However, CRN 
employed thin-window multiwire proportional chambers to detect transition radiation, and 
hence, needed a heavy pressurized container. The proportional tubes used in TRACER 
instead can easily withstand external low pressure; and thus, a pressurized gondola is not 
required and not used for TRACER. 
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3. Balloon Flight and Instrument Performance 

TRACER had first a one-day balloon flig ht in Fort Sumner, Ne w Mexico, in 1999. The 



results from this flight have been published (IGahbauer et al.l 120041 ) . A circumpolar long- 



duration balloon (LDB) flight was originally planned for a flight in the northern hemisphere 
with a launch from Alaska, but could not be conducted because of the lack of relevant in- 
ternational agreements. Therefore, an LDB flight had to be planned in Antarctica, but had 
to wait until 2003, when a new launch vehicle commensurate with the weight of TRACER 
(6,000 lbs., including ballast, telemetry, and balloon-related instrumentation) became oper- 
ational. TRACER was launched near McMurdo, Antarctica, on December 12, 2003. During 
its two-week flight (see figure [3]) at an average height of 37.75 km (3.9 g/cm 2 ) with excur- 
sions not exceeding ±1.75 km, it collected data from a total of 5 x 10 7 cosmic-ray events. 
The data were transmitted at 500 Mbit / sec to the ground (when the instrument was within 
telemetry range) and stored on six disks on-board and recovered after termination of the 
flight. 

Overall, the instrument performed very well during this flight. In particular, the entirely 
passive thermal insulation kept the temperatures of all detector components close to room 
temperature, with diurnal variations of at most a few degrees C. There was no indication 
of a deterioration of the performance of the proportional tubes due to gas poisoning for the 
duration of the flight. This result would be of great practical importance if flights of such 
devices for much longer duration (for instance in space) are anticipated: there seems to be 
no need for gas-purging at the time scales of at least a few weeks, except for the provision 
of make-up gas to correct for minor gas leaks. The data acquisition system operated very 
efficiently, with a dead time of just 6%, but the total time during which data could be 
recorded was limited to ten days, due to a catastrophic failure of the rechargeable lithium 
batteries (which were used to buffer the solar-power system) towards the end of the first 
orbit of the balloon around the South Pole. 



4. Data Analysis 

4.1. Trajectory Reconstruction 

The first step in the analysis of the data is the accurate reconstruction of the trajectory 
of each cosmic-ray particle through the instrument. The entire array of proportional tubes 
is used for this purpose. It consists of four double layers of tubes in each the x- and y- 
direction. Hence, ideally a trajectory is characterized by eight tube "hits" in each projection. 
In practice, sometimes the particle may pass between adjacent tubes of a given layer and 
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not generate a hit in this layer. In addition, spurious signals in tubes outside the particle 
trajectory may be recorded. 

As a first estimate, the path is obtained using the center of each of the tubes hit in an 
event. All possible combinations of these tubes are fit with straight lines in the X and Y 
projection, and the combination in each projection with minimum \ 2 P er degree of freedom 
and maximum number of tubes used is kept. This method allows an accuracy of 5 mm in 
track position. The procedure is quite efficient: less than 5 % of cosmic-ray nuclei are missed 
in the process. 

A preliminary estimate of the charge of each particle is then performed using the tech- 
nique described in the following section. A subset of the data from which we expect the same 
average signal in any tube (high energy O or iVe nuclei) is then used to correct for deviations 
of the tube positions from their "ideal" positions, and to match signal gains within 5 % or 
less. Corrections due to variations in the tube gain with time either due to temperature 
changes or small gas leakages are also calculated and applied. In general, gain variations in a 
given tube are negligible along that tube. Finally, the tube signals are corrected to account 
for minor non-linearities in the front end electronics at large signals. All these corrections 
are necessary to reduce the systematic uncertainties in the tube signals below the level of 
unavoidable physical fluctuations ( "Landau fluctuations" ) in the energy deposited in the gas. 

As a final step, the trajectory is refined by using the fact that the energy deposit in 
each tube is proportional to the track length in that tube. Taking this fact into account, 
an accuracy of 2 mm in the lateral track position is achieved, which corresponds to a 3 % 
uncertainty in the total pathlength through all th e tubes. The accuracy a nd efficiency of this 



method are verified with a GEANT4 simulation (lAgostinelli et al.l 120031 ) of the instrument. 



The accurate knowledge of the particle trajectories is essential to the analysis for two 
reasons: first, it permits corrections of the scintillator and Cherenkov signals due to spatial 
non- uniformities and zenith-angle variations (see section |4~2]) . Second, it makes possible to 
determine accurately the average energy deposit of a cosmic-ray particle when it traverses 
the tubes. Specifically, the average signals are as follows: 

In the top four double layers, the average specific ionization is defined as 



,dE JXlAEj 
dx> YZlAxi 



<-> = tt^ (1) 



and in the bottom four double layers (which form the TRD), the average signal is 



9 



ax 



e;=9 6 a^ 



(2) 



Here, AEi and Axi refer to energy deposit and pathlength, respectively, in tube number 
i. The two quantities defined in (OQ) and (j2j) are, of course, identical within fluctuations at 
energies where no transition radiation is detectable (TR = 0). In order to exclude large 
fluctuations associated with short pathlengths, we remove tubes with Axi < 1 cm from the 
summation. 



As a first step to determine the charge of each particle, the signal recorded by each PMT 
is corrected to account for the spatial non-uniformities in the counter responses. These cor- 
rections are derived using the tracking information together with response maps recorded 
with muons at ground and verified with the flight data themselves in an iterative procedure. 
Subsequently, for each event the average signal of all PMTs in each counter is calculated, 
and then normalized to vertical incidence. Figure H] shows a scatter plot of the top scin- 
tillator signal vs the Cherenkov signal. It is apparent that particles with the same charge 
are clustered along lines. The position along the line depends on the primary energy: a 
concentration of events around maximum Cherenkov and minimum scintillation signals cor- 
responds to energies of minimum ionizing particles and above, and the energies decrease for 
events towards the left of this concentration. 

The scintillator signals depend on particle charge as Z im , due to non-linear effects in the 
light yield, while the Cherenkov signals follow the expected Z 2 dependence. The residual 
scintillation in the Cherenkov counter is visible in figure HI and is subtracted out in the 
analysis. Figure [5] shows a charge histogram for all charges obtained from summing along 
lines of constant charge in figure HI The charge resolution evolves from 0.25 charge units 
for Z=8 (O nuclei) to 0.5 charge units for Z=26 (Fe nuclei) at low energies (< 10 GeV per 
amu), and increases slightly to 0.3 and 0.6 charge units, respectively, at the highest energies. 

Note that figure [5] is generated just on the basis of the signals of the top scintillator, 
and thus depicts a histogram of Z top . One can, in the same fashion, construct a histogram 
of Zb ot) using the signals of the bottom scintillator. For the further analysis, we require that 
Z top and Z bot are consistent within fluctuations for each accepted event. In practice, this 
means that particles with Z bot < Z top are rejected. In almost all cases these are particles 



4.2. Charge Analysis 
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that have undergone a nuclear spallation while traversing the detector material between top 
and bottom of the instrument (for numerical detail, see section 5 and figure [TUl) . 

The procedure to select a sample of a given element for further analysis is as follows: first, 
with a very tight cut on Z botl clean, penetrating nuclei are selected, and their distribution 
in Z top is investigated. This yields the values for charge resolution mentioned above, and 
it indicates the efficiency if a charge cut on Z top is made. If the cut on Z top is of the order 
of AZ top > ±0.5 charge units, the efficiency is quite high. For certain elements, a tighter 
cut on Z top must be made, and the efficiency is then reduced. Numerical details are given 
in section 5 and table [TJ It should be noted that the tight cut on Z bot is only applied in 
order to determine the form of the charge distribution on Z top . For the final analysis, only a 
very loose cut on Z bot is used to remove interacting particles. This cut leads to a negligible 
reduction in overall charge selection efficiency 

Finally, it should be noted that "edge" effects can compromise the charge resolution. 
To avoid this problem, tracks with an impact point within 1 cm of the wavelength-shifter 
bars are excluded. This creates a dead area in the detector of about 4 % of the total area. 

4.3. Energy Analysis 

4-3.1. Overview 

The energy of each cosmic-ray nucleus is obtained from the combined signals of the 
Cherenkov counter and of the proportional tubes. The top four double layers of tubes 
measure the specific ionization of each particle, (dE/dx), while the bottom four double 
layers detect transition radiation, measuring (dE/dx + TR). While the main objective 
of the Cherenkov counter is to identify low energy particles, below the level of minimum 
ionization, the rapid increase of the signal with energy provides a good energy measurement 
up to about 3 GeV per amu. 

At energies between minimum ionization and the onset of TR (3 to 400 GeV per amu) the 
signals in the dE/dx tubes and the TR tubes are the same and increase logarithmically with 
energy. This increase is slow and can be used for an energy measurement unless it is obscured 
by statistical signal fluctuations. The relative level of the fluctuations decreases with the 
nuclear charge as \jZ. Therefore, the energy measurement improves with increasing Z. 
Note that each individual proportional tube provides an estimate for the specific ionization 
by measuring the ratio A^/Axj. Hence, for any particle up to 16 tubes provide independent 
measurements, which are required to be consistent within fluctuations. For the final analysis, 
the average (dE/dx), as defined in equation [1] and [2], is used. 
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At the highest energies (> 400 GeV per amu), the signals from the dE/dx tubes and 
from the TR tubes diverge, and the rapid increase of the TR signals with particle energy 
provides an excellent energy measurement. However, particles in this energy region are 
extremely rare, less abundant than particles in the minimum-ionization region by more 
than four orders of magnitude. To uniquely identify these particles, it is necessary but not 
sufficient to just require that low-energy nuclei are rejected on the basis of their Cherenkov 
signals. In addition, we require that the measurement of their ionization energy loss (dE/dx) 
(performed with four double layers of dE/dx tubes) places them at an energy level well above 
minimum ionization (for details, see section 14.3.31 and table [2]). Thus, the combination of 
dE/dx tubes with the TRD tubes is crucial for the success of the TRACER measurement at 
the highest energies. 



4-3.2. Identification of Sub-Relativistic Particles (< 3 GeV/amu) 

To identify low energy particles, the signals from the Cherenkov counter are compared 
with those of the dE/dx tubes. The dE/dx response, shown in figure [U, is well described by 
the Bethe-Bloch formula. The Cherenkov counter response (figure [1]) must be slightly mod- 
ified to account for effects of delta rays generated by the primary particle in the instrument, 
while the signals from the proportional tube array remain unaffected by d elta rays. These 



effects are well understoo d and have been previously studied and reported (IGahbauer et al. 



20031 ; iRomero- Wolil 120051 ) . The scatter plot in figure [6] shows the correlation of the ioniza- 
tion and Cherenkov signals measured for iron nuclei. The solid line represents the average 
response expected from simulations. The Cherenkov signal increases with energy until it 
reaches saturation, while the ionization signal clearly exhibits the level of minimum ioniza- 
tion (~ 3 GeV per amu). Note that the gray scale of the scatter plot is logarithmic. Plots 
like this one for each charge exhibit the level of fluctuations in the Cherenkov signal: it is 
about 8 % for O nuclei and 2.5 % for Fe nuclei. A cut in the Cherenkov signal is placed at 
a level that corresponds to the minimum in the ionization signal. Events with Cherenkov 
signals below this cut are used to generate their energy spectrum from 0.8 to 2.3 GeV per 
amu. Events with Cherenkov signals above the cut are used to construct the energy spectra 
above 12 GeV per amu as described in the following section. 



4-3.3. High Energy Particles 



After removing all particles below about 3 GeV per amu, the rare high energy parti- 
cles are identified from the combination of the signals from the dE/dx tubes and the TRD. 
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These detectors have been calibrated at accelerators using singly charge d particles. Fig- 
ure [TJ shows the measured response for the TRD of the CRN instrument (jL'Heureux et al" 



19901 ). The radiator combination and gas mixtures of TRACER are identical to those of the 
CRN instrument. However, CRN used plane multiwire proportional chambers rather than 
the layers of cylindrical proportional tubes of TRACER. We have ascertained in Monte- 
Carlo simulations that the accelerator calibrations for CRN remain valid for the detector 
geometry of TRACER . Most of the data shown in figure [7] have been published previously 



(IL'Heureux et al.lll990f ). 



Below 400 GeV/amu we expect no observable transition radiation. Both signals, (dE/dx) 
and (dE/dx + TR), increase logarithmically with energy and will cluster along the diagonal 
in a correlation plot of TRD response vs ionization response as illustrated in figure [SJ Above 
400 GeV per amu, TR becomes observable, and the signal from the TRD will increase above 
the ionization signal. This manifests itself as a deviation from the diagonal in the correlation 
plot. 

As an example, figure M shows the observed cross correlation between the TRD and 
the ionization signals for neon nuclei (Z = 10) above minimum ionization. The small black 
points represent the numerous events with energies below the onset of TR. The rare high 
energy particles with clear TR signals are highlighted. As expected, the data follow the 
response illustrated in figure Note that the highest energy events (the TR events) stand 
out without any background in other regions of the scatter plot. The most energetic neon 
nucleus in this sample of data has an energy of 6 x 10 14 eV. 

We must emphasize again that the off-diagonal position of these higher energy events 
defines them uniquely as TR-events for the selected charge (Z=10 in figure [9]). Any spillover 
due to misidentified charges, with Z either larger or smaller than the selected value, would 
lead to signals along the diagonal, but would not contaminate the well defined "TR-tails" . 



5. Absolute Flux Measurements 

Each event that passes the data analysis cuts is classified in energy as either as a 
Cherenkov Event (< 3 GeV per amu), dE/dx Event (10 - 400 GeV per amu) or TR Event 
(> 400 GeV per amu). Events are sorted into energy bins of width AEi and a differential 
energy spectrum is constructed for each elemental species. We present the spectra in terms 
of an absolute flux dN/dE at the top of the atmosphere. To convert from the number of 
events A A, in a particular energy bin AEi to an absolute differential flux dN/dE(i) one 
must compute the exposure factor, effective aperture, efficiency of the cuts, and unfold the 



13 



instrument response: 




ANi 1 1 1 



(3) 



AE t ' T, ' Si' A 



with ANi the number of events, AEi the energy range, T\ the live-time, e$ the efficiency 



of analysis cuts, A4 the effective aperture and C{ the "overlap correction" due to misidentified 
events from neighbouring energy bins. For Cherenkov events one also must take into account 
that AEi, the energy interval on top of the atmosphere, is not equal to the measured energy 
interval because of energy losses in the atmosphere and in the detector. 

Livetime — The live time is determined by comparing the raw trigger rate with the 
number of events stored on disk. With a raw trigger rate of 60 Hz, the dead time was 
negligible throughout the flight, and the total live time is 787,200 seconds. 

Tracking Efficiency - - The efficiency of the trackfit is 95 % for all charges, and is 
independent of energy. Included is a consistency check in the signal of the tubes contributing 
to the (dE/dx) and (dE/dx + TR) measurements. 

Charge Efficiency — For the more abundant elements (O, Ne, Mg, Si) a charge cut of 
Z top ± 0.5 is made. The efficiencies range from 85 % to 72 % (O, Si) and are listed in table [U 
For the less abundant elements (S, Ar, Ca) a more conservative cut of Z top ± 0.1 is imposed 
to avoid contamination from adjacent charges, which yields a reduced efficiency of 17 %. For 
iron, a cut of Z top ± 1.0 is used at an efficiency of 90 %. The consistency cut on the charge at 
the bottom of the instrument Zb ot can be sufficiently loose to not affect the overall efficiency. 

Effective Aperture — The overall geometric factor of TRACER is 5.04 m 2 sr. For the 
flux computation an effective aperture A4 is calculated, which includes the geometric factor, 
losses due to nuclear spallations, and "dead" regions in the detector. 



with A the area of the detector (206 cm x 206 cm), Pi (9) the probability of survival 
against spallation in the atmosphere and in the instrument as a function of zenith angle 9 (see 
figure [TO!) . and Pd{9) the probability that a particle passing through the top scintillator will 
traverse the entire instrument without encountering any insensitive detector regions. The 
quantity Pd(9) is determined from a Monte Carlo simulation. Insensitive detector regions 
lead to a reduction in the overall geometric factor by about 10 %. 



Ai = A • 2tt 



/ P / (0)P D (0)cos0 d(cos0) 



(4) 



J 8=0 



Losses due to nuclear interactions are computed using a combination of measured 
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cha rge- changing cross se ctions and parameterizations obtained by iHeckmann et al.l (119781 ) 



and IWestfall et al.l (I1979I ). The interaction losses are taken to be independent of energy but 
do depend on the particle mass. Interactions in the atmosphere are determined by measuring 
the residual atmosphere above the payload, which for this flight was monitored by pressure 
sensors and GPS modules attached to the instrument. The average residual atmosphere was 
3.91 g/cm 2 . Interaction probabilities in the instrument are obtained by numerical compu- 
tation using a complete list of materials in TRACER. The fraction of interacting nuclei in 
the instrument can also be obtained experimentally from a comparison of the reconstructed 
nuclear charges for each element in the top and bottom scintillators. Both methods lead to 
consistent results. As an example, an oxygen or iron nucleus at an incident angle of 30° has 
a probability of survival in the atmosphere of 83 % or 73% and probability of survival in the 
instrument of 70 % or 49 %, respectively (see figure [TU1 for all elements). 

Energy Overlap Corrections — As TRACER has three detectors, each with a different 
energy response and energy resolution (see figure [1] and figure [11]) , the widths of the energy 
intervals into which the data are sorted must be appropriately chosen. The intrinsic relative 
signal fluctuations decrease with increasing Z, essentially as 1/Z (figure fTTl) . Therefore, the 
number of acceptable energy bins increases for elements with higher charge number. 

As the shape of the energy spectra is not expected to vary rapidly with energy, choosing 
too fine a subdivision of the energy scale might just lead to increased statistical uncertainty 
in each energy bin, and hence, not improve the quality of the results. More importantly, 
the widths of the bins must be commensurate with the energy resolution of the detector. 
Otherwise, overlap corrections for events being misinterpreted from a neighbouring energy 
bin would become substantial and could possibly bias the result for a steeply falling spectrum. 

In the current analysis of the TRACER data, the bin widths are chosen fairly wide, such 
that the overlap correction factor can be determined with the help of an iterative Monte- 
Carlo simulation. The corrections are computed for each energy bin and element. They are 
fairly small, in general < 20 % (i.e. 0.8 < d < 1.2). 



6. Measured Energy Spectra 

The analysis yields, for each species, numbers of events AAj sorted into energy bins 
AEi = E i+ i — Ei. The ratio AA^/AE'j defines the differential intensity dN/dE for this 
interval: 

dN AN t 

dE {t) = AK (5) 
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We assume that the energy spectrum is represented by a power law CE a over this 
interval, such that: 



AN { = / CE~ a dE 
We then plot the spectral value at an energy level E defined such that: 



(6) 



(7) 



E 



E, 



'I -a) 



(Ell? - Et°)- Xfa 



(8) 



Since E depends on the value of a, an iterative procedure has been used. It is found 
that the dependence of E on a is not very strong. 

For the spectra presented in this work, the highest energy intervals are integral bins, 
with E i+ i — > oo. We then plot the data point at the median energy for this bin, E int : 



E mt = [E]-«/2] 



il/1-a 



(9) 



The value of the flux at this energy is: 



dN 



dE 



AN A a - I) 



int 



E 



X-a 



■E 



int 



(10) 



The energy spectra, in terms of absolute intensities, for the elements O, Ne, Mg, Si, 
S, Ar, Ca, and Fe are presented in tables [2] and [3] and are plotted in figure [12j The low 
energy data points (< 2.3 GeV per amu) come from an analysis of the Cherenkov counter, 
using a restricted angular acceptance. Note that the energy in the table is given in units of 
kinetic energy per amu, whereas in figure [12] the units are given in terms of kinetic energy 
per particle. Data from the TRACER 2003 flight are indicated by the solid squares. For 
clarity the intensity of each element is scal ed by a factor shown o n the left. Existing data 



from measurements in space with HEAO-3 (jEngelmann et al.lll990l ) and CRN (IMuller et al. 



199ll ) are shown for comparison. In the energy regions where overlap exists, the agreement 
of the results appears to be quite good. 
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We note the large range in intensity (ten decades) and particle energy (four decades) 
covered by TRACER. As can be seen, the energy spectra extend up to and beyond 10 14 
eV. The eight elements for which spectra are shown in figure [12] represent the major pri- 
mary species from oxygen to iron; contributions to these species from spallation produced 
secondary elements are expected to be small. We may fit the high-energy part (> 20 GeV 
per amu) of each spectrum to a power law in energy E~ a . The resulting spectral indices 
are shown in figure [13] Remarkably, there is no significant trend of the spectral indices with 
charge Z, and all indices fit well to an average of a = 2.65 ± 0.05. It must be emphasized that 
this result is not an artifact of the method of defining the differential intensities described in 
the previous sections. The similarity of the spectral shapes must indicate that acceleration 
and propagation mechanisms are essentially the same for all species. However, the statis- 
tical errors do not exclude the possibility of deviations from a pure power-law behavior for 
individual spectra. 



7. Comparison with other data 



In figure [12], the energy spectra obtained with TRACER are shown together with mea- 
surements in space with the HEAO-3 satellite up to ~ 40 GeV per amu, and with the CRN 
detector on the Space Shuttle up to about 1000 GeV per amu. However, a multi-decade 
plot such as that of figure d2] may obscure the finer details in a steeply falling spectrum. To 
illustrate this, we show the energy spectra for O, Ne, and Fe, in figure [TH multiplied with 
E 2 - 5 . It is apparent that the differences between individual data sets seem to be mostly 
statistical. 



We note that figure also includes a recent determination of the high-energy iron spec- 
trum that was performed with the HESS Imaging Air Cherenkov Telescope on the ground. 
Here it was possible for the first time to observ e directly the Cherenko v light of iron-group 
nuclei in the atmosphere before they interact (lAharonian et al.ll2007al ). Again, the agree- 
ment with the TRACER data is good, although uncertainty about the choice of a nuclear 
interaction model affects the interpretation of the HESS measurement. 

There are also several recent or current balloon-borne detectors that are designed to 
obtain measurement s of the high-energy cosm ic-ray composit ion, for instan ce RUNJOB 
JPerbina et alil2005h . ATIC JPanov et al-lboOfih . and CREAM JSeo et al.ll20o3 ); all of these 
employ mostly calorimetric techniques. The RUNJOB experiment has published energy 
spectra for mostly elemental groups. The iron spectrum from this experiment is also in- 
cluded in figure [HI Again we see that there is good agreement with the TRACER data. 
Spectra for a few individual heavy elements, more limited in energy coverage and statis- 
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2006; 


Ahn et al. 


2007; 


Zei et al. 


2007) 



labeled preliminary. However, the agreement in the regions of overlap appears to be quite 
good. 

Data for oxygen and iron from TRACER are compared in figure [15] wi th spectra de 
rived from indirect air-shower observations by the EAS-TOP collaboration (INavarra et al. 
200 31), and by the KASCADE group for two different nucleus-nucleus interaction models 
( lAntoni et ajj [2005[ ) . These groups do not report results for individual elements: the fluxes 
for the "CNO group" most likely have about twice the intensity than oxygen alone, while 
the "iron group" probably is dominated by iron. Our results do not yet overlap with the 
energy region of the air shower data, but the gap is becoming smaller, in particular for 
oxygen. Additional measurements should indeed lead to significant constraints on the air 
shower interpretations. 



8. Conclusions 

The energy spectra of the cosmic-ray nuclei determined in this investigation extend to 
energies around 10 14 eV per particle, or significantly higher for the more abundant nuclei. We 
believe that this data set represents the most comprehensive information about the spectra 
of the heavier primary nuclei currently available. 

We have studied in some detail how these results provide constraints on our current 
understanding of the origin of galactic cosmic rays and of models of acceleration and propa- 
gation. We realize that the upper energy limit of the measurement is still below the region of 
the cosmic-ray "knee" , so perhaps it is not surprising that no significant changes in cosmic- 
ray composition are yet observed. However, the great similarity between the individual 
energy spectra is remarkable. The result of our studies of these features is currently being 
prepared for a separate publication. 

It is clear that an extension of the measurements to still higher energy is desirable. The 
dynamic range of the TRD used here is not exhausted with the current measurement. Hence, 
in order to reach higher energy, the TRACER detector would either have to be increased in 
size, which seems impractical as the instrument is already the largest balloon-borne cosmic- 
ray detector in existence, or would have to be subjected to additional long-duration balloon 
flights. Ideally, an instrument like TRACER should be exposed in space for several years. 

An important objective for measurements beyond the current TRACER results would 
be an extension of the dynamic range of the detector, such that the light secondary elements 
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(Li, Be, B) as well as carbon and nitrogen can be covered together with the heavier nuclei. 
If successful to sufficiently high energies (~ 10 3 GeV per amu), this measurement would lead 
to a determination of the L/M abundance ratio, i.e., to a determination of the propagation 
pathlength at high energy. 

This would be an essential step towards a better understanding of the galactic propa- 
gation of cosmic rays. To make this measurement possible, the electronics of TRACER has 
recently been completely refurbished, and another successful long-duration balloon flight was 
performed in 2006. The analysis of these new data is currently in progress and results will 
be reported in due course. 
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Table 1: Charge selection efficiencies 

Element 8 10 12 14 16 18 20 26 

Cut (Charge Units) ±0.5 ±0.5 ±0.5 ±0.5 ±0.1 ±0.1 ±0.1 ±1.0 

Efficiency 0.89 0.83 0.75 0.72 0.17 0.18 0.18 0.90 
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Table 2: Measured Intensities (Oxygen - Silicon) 



Element Energy Range Kinetic Energy 1 Number of Differential Intensity 

(GeV/amu) E (GeV/amu) Events (m 2 s sr GeV/amu)- 1 



O (Z = 8) 



0.8 - 


1.0 


0.9 


25437 


(1.71 


± 0.01) x 


10° 


1.0 - 


1.3 


1.1 


24237 


(1.38 


± 0.01)x 


10° 


1.3 - 


1.9 


1.5 


38226 


(8.95 


± 0.05)x 


HT 1 


12 - 


292 


48 


456918 


(7.27 


± 0.01)x 


io~ 4 


1315 - 


2627 


1834 


77 


(4.4 


± 0.5) x 


io- 8 


2627 - 


6371 


4001 


29 


(5.8 


± 1.1) X 


io- 9 


6371 - 


9834 


7872 


5 


(9.1 


+6. Ox 
-3.9) x 


io- 10 


9834 - 




15611 


6 


(2.0 


+ 1.2N 
-0.8) x 


io- 10 



0.8 - 


1.0 


0.9 


1.0 - 


1.3 


1.1 


1.3 - 


1.9 


1.5 


12 - 


292 


48 


895 - 


1316 


1080 


1 31 6 - 

lOlU 


2627 


1 834 


2627 - 


6371 


4001 


6371 - 


9834 


7872 


9834 - 




15611 


0.8 - 


1.0 


0.9 


1.0 - 


1.3 


1.1 


1.3 - 


1.9 


1.5 


12 - 


292 


48 


696 - 


1316 


946 


1316 - 


2627 


1834 


2627 - 


6371 


4001 


6371 - 




10113 


0.8 - 


1.0 


0.9 


1.0 - 


1.3 


1.1 


1.3 - 


1.9 


1.5 


12 - 


73 


28 


73 - 


292 


138 


567 - 


696 


627 


696 - 


1316 


946 


1316 - 




2088 



(3.08 


± 0.05) x 


10" 


-1 


(2.42 


± 0.04) x 


10" 


-1 


(1.48 


± 0.02) x 


10" 


-1 


(1.21 


± 0.01) x 


10" 


-4 


(3.2 


± 0.8) x 


10" 


-8 


(8.2 


± 2.3) x 


10" 


-9 


(8.9 


+6.9\ 
-4.2) x 


10" 


-10 


(2.2 


±S:g) >< 


10" 


-10 


(7.8 




10" 


-li 



Ne (Z = 10) 0.8 - 1.0 0.9 4031 

3790 
5639 
57991 
16 
13 
4 
1 
2 

Mg (Z = 12) 0.8 - 1.0 0.9 4603 (4.02 ± 0.06) x 10" 1 

4184 (2.96 ± 0.05) x 10" 1 

6339 (1.89 ± 0.02)x 10" 1 

59374 (1.60 ± 0.01) x 10" 4 

23 (3.6 ± 0.7) x IO" 8 
9 (6.3 ± 2 ;?) x IO" 9 
1 (2.4 tli) x IO" 10 
1 (6.7 + 15 5 4 ) x 10- 11 

Si (Z = 14) 0.8 - 1.0 0.9 3127 (3.07 ± 0.06)x IO" 1 

2812 (2.27 ± 0.04)x IO" 1 

4292 (1.39 ± 0.02)x 10" 1 

36369 (5.23 ± 0.03) x IO" 4 

2278 (7.2 ± 0.2) x 10~ 6 

18 (1.3 ± 0.3) x IO" 7 

24 (4.0 ± 0.8) x 10~ 8 
6 (2.2 tH) x IO- 9 



1 see text for definition of E 
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Table 3: Measured Intensities (Sulphur - Iron) 



Element 



Energy Range 
(GeV/amu) 



Kinetic Energy 1 Number of 
E (GeV/amu) Events 



Differential Intensity 
(m 2 s sr GeV/amu) 1 



S (Z = 16) 



Ar (Z = 18) 



0.8 
1.0 
1.3 
1.8 
12 
41 
130 
529 
792 



1.0 

1.3 

1.8 

2.3 

41 

130 

412 

792 



0.9 
1.1 
1.5 
2.0 
21 
70 
223 
644 
1257 



156 
158 
170 
106 
1212 
212 
32 
4 
3 



0.8 - 1.0 


0.9 


70 


(4.0 ± 0.5) x lO- 2 


1.0 - 1.3 


1.1 


61 


(2.9 ± 0.4) x 10~ 2 


1.3 - 1.8 


1.5 


66 


(1.5 ± 0.2) x 10~ 2 


1.8 - 2.3 


2.0 


39 


(8.9 ± 1.4) x 10~ 3 


12 - 41 


21 


454 


(7.1 ± 0.3) x 10~ 5 


41 - 130 


70 


85 


(3.8 ± 0.4) x 10~ 6 


130 - 412 


223 


11 


(1.5 ± 0.5) x 10~ 7 


452 - 696 


558 


2 


(1.9 + H) x 10- 8 


696 - 


1105 


1 


(3.2 tU) xlO- 9 



7.5 ± 0.6) x 10~ 2 
6.3 ± 0.5) x 10~ 2 
3.3 ± 0.3) x 10~ 2 

2.0 ± 0.2) x 10~ 2 

2.1 ± 0.1) x 10~ 4 

1.1 ± 0.1) x 10~ 5 

5.2 ± 0.9) x 10~ 7 
2.1 t\i) x 10~ 8 

2.6 til) xlO- 9 



Ca (Z = 20) 



0.8 - 1.0 


0.9 


99 


(5.5 ± 0.6) x 10~ 2 


1.0 - 1.3 


1.1 


79 


(3.7 ± 0.4) x 10~ 2 


1.3 - 1.8 


1.5 


99 


(2.2 ± 0.2) x 10~ 2 


1.8 - 2.3 


2.0 


61 


(1.3 ± 0.2) x 10~ 2 


12 - 41 


21 


610 


(9.9 ± 0.4) x 10~ 5 


41 - 130 


70 


88 


(4.1 ± 0.4) x 10~ 6 


130 - 412 


223 


19 


(2.8 ± 0.6) x 10~ 7 


452 - 696 


558 


3 


(1.8 tl'l) x 10- 8 


696 - 


1105 


1 


(1.9 ±f"l) x 10- 9 



Fe (Z = 26) 



0.8 
1.0 
1.3 
1.8 
12 
23 
41 
92 
529 
606 
895 



1.0 
1.3 
1.8 
2.3 

23 

41 

92 

412 

606 

895 



0.9 
1.1 
1.5 
2.0 
16 
30 
60 
183 
566 
733 
1421 



1821 
1669 
2166 
1388 
10589 
4468 
2132 
543 
9 
12 
11 



(2.21 
(1.65 
(1.02 
(6.2 
(1.22 
(2.80 
(4.7 
(1.9 
(1.3 
(5.2 
(7.4 



± 0.05) x 
± 0.04) x 
± 0.02) x 
± 0.2) x 
± 0.02)x 
± 0.04) x 
± 0.1) x 
± 0.1) x 

+0.6N 
-OA) x 

± 1.5) x 
± 2.2) x 



io- 1 
io- 1 
io- 1 
io- 2 

IO" 3 

io~ 4 

IO" 5 
IO" 6 

io- 7 
io- 8 
io- 9 



1 sce text for definition of E 
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Fig. 1. — Energy response of detection techniques used in TRACER: Cherenkov counter 
(CER), energy loss in gases (DEDX), and Transition Radiation (TRD). 
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Fig. 2. — Schematic drawing of TRACER 
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Fig. 3. — Trajectory of 2003 TRACER Antarctic Long Duration Balloon Flight. 
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Fig. 4. — Scatter plot of top scintillator vs. Cherenkov signals in arbitrary units. 
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Fig. 5. — Charge histogram for all events measured in flight. 
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Fig. 6. — Scatter plot of dE/dx vs. Cherenkov signals for iron nuclei. The black line is the 
average response obtained from simulations. The gray scale is logarithmic. 
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Fig. 7. — Energy response of the Transition Radiation Detector as measured at accelerators. 
Excep t for two highest-energ y data points, these measurements have been published previ- 
ously (IL'Heureux et al.lll990l ). The open symbols and dashed lines refer to measurements of 
dE/dx only. 
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SQRT (Specif ic Ionization) 

Fig. 8. — Correlation of average responses (square root) of Transition Radiation and specific 
ionization detectors for relativistic nuclei. Four elements are displayed to illustrate the charge 
dependence of the responses. 
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Fig. 9. — Scatter plot of TR vs. dE/dx signal for neon nuclei. The highlighted points 
represent the highest energy events measured with the TRD. As expected the transition 
radiation events have signals in the dE/dx detector which are well above the minimum 
ionization level. Units are arbitrary. 
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Fig. 10. — Survival probability all nuclei incident at a zenith angle of 30° in a residual 
atmosphere of 3.91 g/cm 2 (dash-dot), within the TRACER instrument (dash), and total 
survival probability (solid). 
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Fig. 11. — Energy resolution (la) of individual detector subsystems vs charge Z, and for 
typical energies: Cherenkov: 7 = 2 (dashed line), dE/dx: 7 = 100 (dash-dot) and TRD: 7 
= 1000 (solid line) 
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Fig. 12. — Differential energy spectra vs energy per particle of the cosmic ray nuclei : O, Ne, 
Mg, Si, S, Ar, Ca and Fe. Results from the TRACER 2003 flight are indi cated by the solid 
square s. Existing data from the HEAO-3 experi ment (open d i amond s : lEngelmann et al. 
(119901 )) and the CRN experiment (open crosses : iMiiller et al.l (119911 )) are shown for com- 
parison. The dashed line represents an independent power-law fit to each spectrum above 
20 GeV per amu. 
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Fig. 13. — Spectral indices of a best power-law fit to the combined TRACER and CRN data 
above 20 GeV per amu. The line indicates the an average spectral fit of E~ 2,65 . 
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Fig. 14. — Energy spectra multiplied with E 2,5 fo r O, Ne, and Fe, for diffe rent observations: 
TRACER (solid sq uares), HEAO-3 ( diamon ds: lEngelmann et all jl990h ). CRN 



crosses: 



Miiller et al.l (Il99lf l). RUN JOB (rilled circles: iDerbina et al.l (1200.1 n and HESS (open stars: 
QGSJET model, filled stars: SYBILL model; lAharonian et al.l (j2007al )). The error bars 
shown are statistical. 
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Fig. 15. — Energy Spectra from TRACER (solid squares) compared with the interpreta- 
tion of air shower d ata of KASCADE (solid triangles, for two different interaction models : 
Antoni et al.l (120051 )) and of EAS-TOP (open cros ses, two data points for each energy repre- 
sent upper and lower limits: iNavarra et al.l (120031 )). The spectra are for oxygen and for iron 
for TRACER, but for the "CNO-group" and the "Fe-group" for the other observations. 



